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Abstract 

In this letter we investigate link scheduling algorithms for throughput maximization in multicast 
wireless networks. According to our system model, each source node transmits to a multicast group 
that resides one hop away. We adopt the physical interference model to reflect the aggregate signal to 
interference and noise ratio (SINR) at each node of the multicast group. We present an ILP formulation 
of the aforementioned problem. The basic feature of the problem formulation is that it decomposes 
the single multicast session into the corresponding point-to-point links. The rationale is that a solution 
algorithm has more flexibility regarding the scheduling options for individual nodes. The extended MILP 
problem that also considers power control is solved with LP relaxation. Performance results for both 
the ILP and MILP problems are obtained for different traffic loads and different number of nodes per 
\f} , multicast group. 
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Index Terms 



Link scheduling, wireless multicast, wireless networks, power allocation, integer linear program, 
^ , mixed integer linear program, approximation algorithm. 

I. Introduction 

Wireless multicasting poses significant technical challenges that have attracted considerable 
amount of research work for several years. The problem is more relevant than ever due to the 
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widespread use of WiFi-enabled mobile devices that need to transmit high volumes of data 
to several users. Smart-phones with multicasting capabilities are envisioned as one of the key 
adopters of such technologies especially for the distribution of high-quality locally captured 
video data (at the multicast source). To this aim there are notable practical efforts for wireless 
multicasting in WiFi networks 0], CD- In this letter we are concerned with an environment 
where multicasting is routinely adopted for wireless transmission among nodes in a distributed 
network (see Fig. [T]). This means that several nodes may need to concurrently use multicast 
communication. 

In wireless multicast one of the key problems is the impact of heterogeneous destinations on 
the performance of the entire multicast session [2]. A destination that belongs to the multicast 
group, but is characterized by the worst channel quality among all the destinations, will be 
the bottleneck of the multicast communication since it will require increased transmit power 
and retransmissions. A scenario beyond a single multicast session needs to consider a wireless 
network where several muticast sessions are active in the same space, time period, or even 
frequency band. Therefore, it is also possible that a multicast session generates interference to 
the destinations of the rest of the multicast sessions. This creates an additional problem for the 
performance of multicast in the wireless network. Fig. [TJ depicts this situation where sources 
Si, S 2 that multicast to their respective groups D n , D 12 , D 13 and D 2 i,D 22 respectively. In this 
case Si must increase the transmit power to a level that the packet is decodable also from Di 2 . 
However, this decision will generate interference to node D 2 i and it may render undecodable the 
second multicast transmission. The same is true for decisions made by S 2 (transmission range 
is not shown). 

The problem we address in the context of the previously described scenario is the following. 
Given a network with a set of next-hop multicast sessions how can we schedule the wireless 
multicast transmissions so that throughput is maximized? Significant performance improvements 
are observed and originate from a re-formulation of the integer linear problem (ILP) that de- 
scribes the throughput maximization problem for this scenario. Subsequently, the ILP problem 
formulation is extended to a mixed ILP (MILP) that also considers power control and is solved 
with LP relaxation and a heuristic that exploits the multicast nature of the problem. The only 
disadvantage of the proposed approach is a higher number of optimization variables. 

Scheduling multicast transmissions in wireless networks has attracted a certain amount of work 
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Fig. 1. Example network with two concurrent multicast sessions. When both nodes multicast during the same slot, then they 
both interfere to at least one of the destinations of the other multicast group (left). When the transmission power of Si is reduced 
and it can only transmit reliably to Dn,Dn then the result is that S% can multicast reliably to both its destinations (right). The 
optimal activations of the links that compose each multicast session is addressed in this work. 

even a few years earlier. Wieselthier et. al had shown in the first results that indicated that 
with proper power allocation (called the broadcast incremental power algorithm) it is possible 
to minimize interference for multicast sessions. Heuristic joint power control and scheduling 
algorithms were also studied by Wang et. al in H|. This work is of particular interest since 
the authors focused on scheduling for wireless multicasting but with the objective to minimize 
power consumption. For a given group of multicast sources and their corresponding destinations, 
the power optimization problem was formulated as a MILP. The optimal values of the transmit 
power were evaluated so that the SINR requirements at the receivers were fulfilled while the total 
power expenditure was minimized. In more recent works Gopinathan et. al presented a model 
for optimal multicast in multi-channel multi-radio wireless networks under the assumption that 
the channel assignment is static 0. Krishnan et. al focused on the problem of identifying 
the optimal multicast trees by considering also the next hops in multi-hop wireless ad hoc 
networks [6J. More theoretical studies driven from an information-theory perspective, focus on 
exploiting the broadcast advantage for multiple sessions 0. In multicasting is considered 
in a system that employs network coding and allows cross-layer interactions while the authors 
follow a utility optimization scheme for calculating the optimal source rate. 

II. System Model and Initial Problem Formulation 

We study a network model where a set S = {S±, S2, SV} of sources want to communicate 
with a number of multicast destination nodes that are denoted as the set V. Each source multicasts 
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to D nodes. The definition of a multicast link in this letter is extended and it involves a number 
of nodes that is equal to the size of the multicast group. We refer explicitly to a link between two 
nodes as a point-to-point link. The complete network is modeled with a directed graph F(U, V), 
where U and V are the set of point-to-point directional links, and the set of nodes, respectively 
(V = {S, £>}). We also use the conflict graph of F that is defined as G(U, V), and contains the 
interfering relationships among the N x D point-to-point links in the network. Each vertex in 
the conflict graph represents a wireless link in the network, and there is an edge between two 
vertices if and only if the links represented by the vertices conflict (i.e. they interfere with each 
other and simultaneous transmission is impossible). On the other hand, a clique in the conflict 
graph represents a group of links that cannot transmit concurrently, and hence they must access 
the channel exclusively. 

A. Initial MILP Formulation 

In the most closely related work to this letter the objective was power minimization for a 
wireless multicast scenario [4]. We follow a similar approach for defining the initial multicast 
optimization problem but in our case we consider throughput maximization. In this formulation, 
the binary optimization variable x\ indicates whether a transmission from the source of the 
multicast group % occurs in slot t. T is the maximum number of slots, Pf is an optimization 
variable that corresponds to the transmit power of source i during slot t, (3 is the SINR decoding 
threshold at the destination, a 2 is the AWGN variance, 7^ = l/d^ where rfy is the distance 
between source i and destination j while a is usually set in a value between 3 and 4. The actual 
problem formulation named MC — ALL is given below: 



j T N 

max — > > Dax ■ 

" 1 t=l i=l 



PhiJ '' n xi)A >p,VjeVi,ieS,teT(i) 



° 2 + Ekes-W P bk 3 

T T 

5>J>l,VieS(2), ^it<ir ai ,VzeS(3), 

t=i t=i 

< P\ < P*' maX 4 (4), **G{0,1}(5) 



The objective is to maximize the throughput by increasing the number of scheduled sources and 
after taking into account as a weight the number of nodes in the multicast group of each source 
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Di. Next, constraint (1) is the SINR constraint. In (1) A is a large number needed for ensuring 
that the SINR constraint is satisfied when the respective link is not scheduled. In our performance 
evaluation it is set to a value higher than the maximum possible SNR in the network. The key 
observation from the formulation in (1) is that it ensures that the SINR constraints are satisfied 
for all the multicast destinations of each source. The remaining constraints (2),(3),(4) ensure 
first that each source is scheduled at least once during the T slots, second that the transmitter 
power is limited according to a maximum value p. ma:E for the complete set of T slots, and third 
that that the transmitter power p^ max for a specific slot is limited depending on the transceiver 
specifications. A note should be made here that will set the stage for our proposed optimization 
approach. Whenever the optimal solution cannot be found for the given set of multicast sources, 
strong interferers are completely eliminated (and this means complete multicast groups cannot 
be scheduled in that slot, i.e. x\ = 0). This is a key disadvantage of the MC — ALL approach. 

III. Throughput Maximization in Multicast Wireless Networks 

A. Proposed MILP Formulation 

The last observation we made for the behavior of MC — ALL stems from one critical detail. In 
the previous formulation a specific multicast group is treated as a single schedulable link/entity. 
Avoiding this limitation is necessary for allowing more flexible scheduling decisions. To this 
aim we relax the requirement that with a single broadcast transmission from a source, all the 
associated multicast destinations must receive the packet. We introduce now the optimization 
variable x\j that indicates the activation of a single point-to-point link from multicast source 
i to a specific destination j. The optimization variables for the transmitter power remain the 
same and are denoted again as P\. The DMC — OPT (Dynamic Multicast) problem allows 
the source to activate dynamically a subset of its point-to-point links that compose a multicast 
link in a specific slot t and not the complete set in order to achieve throughput optimality. 
For the remaining non- scheduled destinations we introduce constraints that ensure that they are 
scheduled during a number of slots. The detailed formulation is named DMC — OPT and is 
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given below: 



T N Di 



• <=i i=i j=i 

^• + (l-4)A 



XX- <A,Vie5(2), 4 g{0,1} (3) 



T 



t=i 

T 

5X < p r x (s), ^' ,mm < ^ < p/ ,ma x (6) 



t=i 



Constraint (1) ensures that each of the D { destinations, that are counted with the subscript j and 
correspond to the multicast group originating from source i, must have the SINR higher than 
the required threshold j3. The next constraint (2) essentially says that when a node multicasts, 
this action corresponds to the activation of a number of source-destination point-to-point links 
that their number must be less or equal to the number of multicast destinations Di. With this 
constraint it is possible that not all point-to-point links are activated. Although a signal from 
a transmission might be received at every network node, with the term activation we mean the 
selection of a proper power level so that the corresponding destination can decode the packet (i.e. 
the SINR is above f3). A valid or 1 value for the activation of a particular link is ensured with 
constraint (3). It is important to clarify constraint (4) that basically ensures that each destination 
of a multicast group, that corresponds to source % is scheduled at least B iy j slots. With these last 
two constraints that we explained it is possible that a multicast group is changing dynamically 
on a slot basis depending on the interference conditions. On the other hand constraint (5) for 
the power budget is also very crucial for ensuring fairness with MC — ALL. This constraint 
ensures that the total power expenditure for a source during the T slots is within a certain budget. 
Therefore, the source must comply with this power budget regardless of how many times it was 
activated for completing a single multicast transmission. Finally, (6) ensures that the transmit 
power constraint per-slot must be between a minimum and a maximum value. 

The MILP is solved by relaxing the integer constraints for the slot activation indicator variables 
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x\a so that a linear program (LP) is solved. The practical problem we have to address is the 
selection from the result of the relaxed LP, the optimal binary result that should either be or 
1. In works such as [8] the randomized rounding method was employed for this purpose while 
similar approaches were followed in [9] . In this letter we exploit the nature of the problem, that 
is the multicast transmission, in order to reach the desired result. The approximation algorithm 
shown in Fig. [2] is described next. 

After the relaxed LP problem is solved (line 1) and the optimal solution is derived in the form 
of the vectors x, P, the algorithm calculates the parameter Cj = J2f=i Ylt=i x tj m at expresses 
the average number of slots that source % must be activated. The main idea of the approximation 
algorithm is to re-order the multicast sources according to the value of q. The source that 
has the highest value for Cj must be scheduled for achieving throughput optimality since its 
transmission can reach the highest number of multicast destinations. To accomplish that, the 
algorithm calculates also which is the maximum number of links that can be active (lines 9- 
10). For the specific source that has the maximum |~Cj] , the corresponding x\^ for each destination 
are set to 1 (lines 1 1-14). If this schedule is feasible the algorithm selects it and moves to the next 
step. Otherwise it takes a number of activated point-to-point links x\^ to be equal to [cjj (line 
10) which is feasible by definition since q corresponds to a feasible solution. The corresponding 
variables x\j are again set to 1 and are stored in the vector x, while the solution for transmit 
power is stored in Pj. The new LP is solved again by using as input the set Z of unscheduled 
sources and the now constant 5q, FV 

B. ILP Formulation with Constant Transmitter Power 

By considering a constant transmitter power from all the sources the problem can be signif- 
icantly simplified leading to considerably interesting results even without power control. The 
problem is defined by removing constraint (6) and by setting P\ to a constant power level P in 
DMC - OPT. The solution to the above problem is obtained with CPLEX 12.04 in the results 
section. 

IV. Performance Evaluation 

In this section we compare the performance of the proposed scheduling algorithms named 
DMC - OPT to that of the multicast scheduling algorithm DMC - ALL where the hyper- 
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milp_relax{S, V, T, B) 
l: x, P = lp(S, V, T, B) //find the solution vectors 
2: for all nodes % G S do 

4: end for 

5: Z = S 

6: for all % e Z do 

7: cat = 

8: while i not feasible do 

9: if {cat == 0) i = argmax(cj), m = |~q] 
10: if {cat == 1) i = argmax(cj), m = [qJ 
11: for p2p link j = 1 until m do 
12: j = arg max 

13: set xlj = l,pP = Pj 

14: end for 

15: if DMC-OPT is feasible then 

16: feasible = TRUE, cat = 

17: else 

18: feasible = FALSE, cat = 1 

19: end if 

20: end while // Solution for i: x;, Pj 

21: Z = Z-{i} 

22: X, P = lp(Z, T>, T, B, Xj, Pj) 

23: end for // The final solution: x, P 

Fig. 2. Pseudocode for the approximation algorithm of the relaxed DMC — OPT MILP problem. 

graph of each mulitcast link is treated as one schedulable entity. We also present results for the 
scheduling of unicast transmissions named UNI — ALL. The SNR threshold j3 is lOdB, the 
maximum and minimum transmit power levels are P-' max =300mW, p^ min = Q$\p^ max > while 
for the ILP case P=0.3P-' max . The node distances d t j are randomly and uniformly selected is 
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Fig. 3. Results for ILP. 



the range [0, 1] with the path loss exponent a=3. Also cr 2 =0.1mW. We consider multicast groups 
with different number of destinations and different traffic loads in terms of the required active 
slots B. 

In Fig. [3] the performance after solving the ILP is presented for both DMC — ALL and 
DMC — OPT with a multicast group of two destinations. In Fig. |3£a) where B=T=S we see 
that the proposed scheme leads to high throughput increase. It is important to observe that 
the performance of the DMC — OPT scheme reaches a peak at a slightly higher number of 
multicast sources/groups. Also note that as the number of sources is increased, the performance 
of DMC — ALL deteriorates faster than the performance of UNI — ALL and this is only 
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Fig. 4. Results for MILP. 

because the increased node density increases interference. For a lighter traffic load of B=A and 
T=8 in Fig. Ob), the performance trend is similar. 

The results for the LP relaxation of the MILP problem are shown in Fig. HI With dashed lines 
we present the optimal solution calculated with CPLEX. We observe in Fig. Ufa) that for a traffic 
requirement of B=4 out of T=8 slots, there is a performance improvement for DMC — OPT 
that is increased as the number of multicast destinations is increased. This behavior occurs 
primarily for a number of sources that is less than the number of maximum number of used 
slots T while for high node density interference dominates again. For a higher traffic load 
and backlogged nodes (£?=T=8) the results can be seen in Fig. Sb). For different number of 
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destinations, the performance improvement of DMC — OPT is higher over DMC — ALL 
even for a small number of sources. However, in both cases and for higher node density the 
performance of all schemes converges. The reason is that interference is higher and fewer options 
for scheduling exist (the activation of a source for one destination generates nearly the same 
interference even if more destinations are activated). From the results that represent the optimal 
solution obtained with CPLEX, we can see that the approximation algorithm is more sensitive 
to the number of multicast destinations than the actual number of sources. For a lighter traffic 
load the proposed algorithm can approach closer the optimal solution. 

V. Conclusions 

In this letter we presented ILP & MILP formulations for the problem of multicast link 
scheduling in wireless networks. With the proposed formulation a single multicast transmission 
is separated across different time slots while complying with the power budget. For the MILP 
formulation we proposed an approximation algorithm that exploits the multicast nature of the 
problem. The performance results indicate that as the multicast group is increased, it is more 
critical to employ the proposed approach that allows the scheduling algorithm to freely allocate 
individual transmissions across time. For constant transmit power, the proposed approach offers 
higher throughput benefits because with existing schemes many multicast links are disqualified 
from being scheduled. In our future work our plan is to derive an analytical approximation bound 
for the proposed algorithm with respect to the optimal solution. 
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